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EFFICIENT EQUALIZATION FOR DETECTION OF SYMBOLS IN DIGITAL 

BURST TRANSMISSIONS 

FIELD OF THE INVENTION 

The invention relates generally to the detection of symbols in digital burst 
transmissions and, more particularly, to equalization of received bursts in order to reduce the 
mean squared error (MSE) of the detected symbols. 

BACKGROUND OF THE INVENTION 

FIGURE 1 diagrammatically illustrates pertinent portions of an exemplary digital 
communication system according to the prior art. In the example of FIGURE 1, a digital 
transmitter produces a sequence of iid (independent identical distributed) symbols designated 
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as A[n] . This sequence of symbols is transmitted in a digital communication burst across an 
equivalent channel to an equalizer 1 1 provided in a digital receiver. As illustrated in 
FIGURE 1, the equivalent channel between the transmitted symbols A[n] and the equalizer 
input includes portions of both the digital transmitter and the digital receiver, as well as a 
physical communication link or channel. For example, the equivalent channel would include 
modulation components of the transmitter, and demodulation components and notch filters 
of the receiver. 

The equalizer 11 includes a feed forward equalizer (FFE) 13, which is a finite 
impulse response (FIR) filter that processes the received samples after some pre-processing 
in earlier stages of the receiver (included in the equivalent channel of FIGURE 1). The 
sampling rate of the input to this feed forward equalizer module can be one sample per 
symbol or more. A sampling rate of one sample per symbol matches exactly the symbol rate 
(also known as the baud rate), and a feed forward equalizer having a one sample per symbol 
sampling rate is conventionally referred to as a T- spaced equalizer (TSE). A feed forward 
equalizer whose input sampling rate is higher than the symbol rate (for example twice the 
symbol rate) is conventionally referred to as a fractionally spaced equalizer (FSE). An 
upsampler 18 illustrates the possibility of multiplying the sampling rate by M at the 
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transmitter, in which case a corresponding downsampler 14 divides the sampling rate by M 
at the receiver. 

A slicer 15 is the decision element of the equalizer. The output of the feed forward 
equalizer module 13 is combined with information from a feedback loop 12 to produce the 
input for the slicer 15. The output Z[n] of the slicer 15 is indicative of which symbol (in the 
constellation) was transmitted. One typical conventional decision criteria that can be used 
by the slicer 1 5 is the Minimum Euclidean Distance (MED) criteria. The output of the slicer 
15 is also provided as an input to the feedback loop 12. 

The feedback loop 12 includes a feedback filter 17, for example a FIR filter that 
operates at the symbol rate. The input to this feedback filter 17 is the slicer output, namely 
the decision regarding what symbol was last transmitted. Assuming correct decisions from 
the slicer 1 5 , the feedback filter 1 7 can compensate for the post cursor part of the interference 
from symbols that are previous in time. Such interference from timewise previous symbols 
is conventionally known as inter-symbol interference or ISI. 

The FIR filter of the feed forward equalizer 13 is defined by a finite length sequence 
C[n] of filter coefficients (or filter taps), and the FIR filter of the feedback filter module 17 
is defined by a finite length sequence D[n] of filter coefficients. The length of C[n] is herein 
referred to as K, and the length of D[n] is herein referred to as L. The basic task of the 
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equalizer is to find the coefficients of C[n] and D[n] that minimize the mean squared error 
(MSE) associated with the decided symbols in the slicer output Z[n]. A coefficient 
determiner 19 uses conventional techniques, for example, least mean squares (LMS), least 
squares (LS) or recursive least squares (RLS) to determine the coefficients of C[n] and D[n] 
that minimize the MSE of the slicer 15. 

FIGURE 2 diagrammatically illustrates a conventional model of the equivalent 
y channel of FIGURE 1 . The equivalent channel is modeled as a digital filter 2 1 having a filter 

Li*: 

tap sequence H[n], and an additive noise factor. Assume for purposes of exemplary 
; S exposition that the additive noise in the equivalent channel model is stationary white noise 

is 

10 ;*£ w[n] that obeys the following statistics: 

□ Equation 1 : 

■J! 

K E{w[n]} = 0 

Equation 2: 

{NO/2 1=0 1 

Rm{1]= E{M{n]-w{n-l]} = \ \ 

15 I 0 otherwise) 
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where E is the expected value operator, NO/2 is the noise variance, R is the autocorrelation 
function, and / is the time index. 
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The source A[n] of FIGURE 1 is also white (i.e., independent identical distributed) 
and obeys the following statistics: 
Equation 3: 

^ v otherwise J 

where rr ^ is the signal variance. 

u A 

Under the exemplary conditions described above, the MSE of the slicer output can 
be expressed as follows for the case of a T- spaced equalizer: 
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For the fractionally spaced equalizer, the MSE can be expressed as follows: 
Equation 5: 
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For a T- spaced equalizer, Equation 4 is the cost function that is to be minimized with 
respect to C[n] and D[n]. In the case of a fractionally spaced equalizer, Equation 5 is to be 
minimized with respect to C[n] and D[n]. 
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In Equations 4 and 5 and hereinafter, U is the number of uncausal coefficients in 
C[n], and V is the number of causal coefficients in C[n]. 

For purposes of exemplary exposition, the following conditions are assumed. 
Communications occur in bursts through the equivalent channel of FIGURE 1, which 
equivalent channel can vary between two consecutive bursts. A single burst includes a 
preamble (a short, unique training sequence) and an unknown sequence of symbols (the 
substantive information), and the channel can be easily estimated from the preamble at the 
beginning of the burst. It is desirable to decode the burst with minimum latency, and the 
number of MMACs (Million Multiply and Accumulate per second) is of course limited. 
Under these conditions, the receiver has a limited time and a limited number of data 
processing operations to achieve a sufficiently low MSE. 

Conventional techniques which can be implemented by the coefficient determiner of 
FIGURE 1 to determine the C[n] and D[n] coefficients which will minimize the MSE 
include the least mean square (LMS) algorithm and its various versions. To use these 
decision directed algorithms, the MSE between the slicer input and the sheer output must 
be lower than a certain threshold (which threshold depends on the constellation being used) 
in order to permit the decoding of the unknown symbols to commence with an acceptably 
low symbol error rate. Most LMS algorithms therefore have the disadvantage of potentially 
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requiring many iterations and many repetitions over the preamble to achieve the desired 
results. 

The conventional LS or RLS algorithm can be used to determine the C[n] and D[n] 
coefficients, but this algorithm is more complex and therefore less practical than the LMS 
algorithm. If the channel is unknown, implementation of the LS/RLS approach involves 
matrix inversion of the size K + L, where K is the length of the sequence C[n] and L is the 
length of the sequence D[n]. If the channel is known, then a matrix inversion of size K is 
required. Although the LS/RLS, algorithm provides the optimal solution (based on the 
empirical statistics of the preamble) in terms of MSE, the aforementioned matrix inversions 
are quite expensive in terms of the required data processing power (e.g. multiply and 
accumulate operations). 

It is desirable in view of the foregoing to provide for equalization that requires less 
computational complexity than the prior art, while still maintaining sufficient performance. 

The present invention provides efficient equalization by initializing the values of only 
a subset of the feed forward and feedback coefficients. After these values are initialized, 
decoding of the unknown symbols is commenced using, for example, the LMS algorithm or 
another iterative, decision directed algorithm to determine all of the feed forward and 
feedback coefficients. The initialization of a subset of the coefficient values before decoding 
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the unknown symbols can advantageously provide an acceptable trade-off between 
computational complexity and equalization performance. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 diagrammatically illustrates pertinent portions of an exemplary digital 
communication system according to the prior art. 

FIGURE 2 diagrammatically illustrates a model of the equivalent channel of FIGURE 

1. 

FIGURE 3 diagrammatically illustrates pertinent portions of an exemplary 
embodiment of a digital communication system according to the present invention. 

FIGURE 4 diagrammatically illustrates pertinent portions of an exemplary 
embodiment of the coefficient initializer of FIGURE 3. 

FIGURE 5 diagrammatically illustrates exemplary operations which can be 
performed by the coefficient determiner embodiments of FIGURES 3 and 4. 



-10- 



102521.00031:2453403.1 



Patent Application 
Docket No. TI-29893 



DETAILED DESCRIPTION 

FIGURE 3 illustrates pertinent portions of an exemplary embodiment of a digital 
communication system according to the present invention. In the embodiment of FIGURE 3, 
the equalizer 30 includes a coefficient determiner 31 having a coefficient initializer 33 and 
a coefficient adaptor 35 and a selector 37 for selectively coupling the coefficient initializer 
and the coefficient adaptor to the feed forward equalizer 1 3 and the feedback filter 1 7. Other 
than the coefficient determiner 31, the components of the exemplary embodiment of 
FIGURE 3 can be generally the same as the corresponding components illustrated in 
FIGURE 1. 

The coefficient initializer 33 determines initial coefficient values of a subset of the 
coefficients of C, and also determines initial coefficient values of the coefficients of D. 
These initial coefficient values are output at 34 by the coefficient initializer 33, and are then 
provided to the feed forward equalizer 13 and the feedback filter 17 via the selector 37. 
Thus, a subset (less than all) of the C and D coefficients is initially provided to the 
equalizer 30. A control signal 36 initially selects the coefficient initializer output 34 to be 
passed through the selector 37, thereby permitting the initial coefficient values to be applied 
to the filters 13 and 17. After the initialized subset of C and D coefficient values has been 
output at 34, the control signal at 36 causes the selector 37 to switch and select the output 38 
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of the coefficient adaptor 35. The control signal 36 can be, for example, a suitable output 
of a state machine (not explicitly shown), which output can become active when a 
predetermined amount of time (adequate for calculation of the initialized subset of C and D 
coefficient values) has elapsed since arrival of a communication burst at the receiver. After 
the control signal 36 becomes active to switch the selector 37, the coefficient adaptor 35 
calculates all or some of the C and D coefficients using, for example, any desired 
conventional iterative, decision directed algorithm such as LMS. The subset of coefficients 
initially provided by the coefficient initializer 33 provides a MSE which is low enough to 
thereafter begin decoding the symbols according to the coefficients calculated by the 
coefficient calculator 35, while still maintaining an acceptable symbol error rate. Thus, the 
coefficient calculator 35 calculates all of the subsequent feed forward coefficients of C and 
the subsequent feedback coefficients of D in generally conventional fashion (analogously to 
the operation of the coefficient determiner 19 of FIGURE 1). The initialized subset of 
coefficient values, i.e. initialized values of a subset of the C coefficients together with 
initialized values of the D coefficients, is provided by the coefficient initializer 33 to the 
coefficient calculator 35 which in turn uses these initial values to begin the calculation of the 
subsequent coefficients of C and D. 
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FIGURE 4 diagrammatically illustrates pertinent portions of an exemplary 
embodiment of the coefficient initializer of FIGURE 3. The embodiment of FIGURE 4 
includes a first coefficient calculator 41 and a second coefficient calculator 43 for 
respectively calculating the initial subset of C coefficient values and the initial D coefficient 
values. The first coefficient calculator 41 calculates the initial subset of C coefficient values 
in response to three inputs, namely an input 42 for receiving channel estimation information, 
an input 44 for receiving linear system information, and an input 45 for receiving weight 
matrix information. In response to the aforementioned information input at 42, 44, and 45, 
the first coefficient calculator 41 implements a weighted least squares solution for a chosen 
subset of the C coefficients. 

The channel estimation information is conventionally available from conventional 
channel estimation procedures performed on the preamble of the communication bursts. 
This channel estimation information corresponds to the filter tap sequence H[n] of 
FIGURE 2. H[n] consists of a sequence of numbers, for example (3, -5, 7, 9, 2, -0.5), which 
numbers are the filter taps of the digital filter component of the equivalent channel. The 
coefficient calculator 41 chooses one of these taps as the leading tap of the sequence H[n]. 
Although any of the numbers in the sequence could be chosen as the leading tap, in this 
example the number having the largest absolute value, namely 9, is chosen as the leading tap. 
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Once the leading tap is chosen, the remainder of the taps of H[n] are indexed as follows: 
uncausal taps are H[-3]=3, H[-2]= -5, and H[-l]=7; leading tap is H[0]=9; and causal taps 
are H[l]=2, and H[2]= -0.5. The aforementioned uncausal and causal taps have the 
following physical interpretations. The uncausal taps mean "physically" that injecting a unit 
symbol (for example numeral 1) to the system (i.e., the channel) at time t=0 will cause an 
output signal at earlier times, namely t= -3, t= -2 and t= -1 . That is, the channel, "responds" 
before it is triggered. The leading tap is the channel output at time t=0, that is, at the time 
that the channel is actually triggered. The causal taps represent the response after the channel 
is triggered. An ideal channel will have only one tap in its response, namely the leading tap, 
and will have no other reflections. The best leading tap is the one that, when chosen together 
with optimal values of the C coefficients and the D coefficients, results in minimum MSE. 

After the leading, causal and uncausal taps of H[n] have been determined, the C 
coefficient calculator 41 finds a weighted least square (WLS) solution for only a selected 
number (less than all) of the coefficients of C. The number N of coefficients to initialize in 
C can be determined, for example, based on factors such as the computational feasibility of 
inverting an NxN matrix, and whether or not N initialized coefficients provides the desired 
performance in terms of initial MSE and symbol error rate. Factors such as these and others 
can be evaluated, for example, in view of empirical information obtained from 
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experimentation and simulation. In one exemplary embodiment, C[n] includes 13 
coefficients, and the C coefficient calculator 41 of FIGURE 4 calculates initial values for 
only 3 of those coefficients. 

In general, the purpose of the N initialized coefficients is to minimize the uncausal 
part of the channel and FFE common impulse response. Any causal part will be handled by 
the feedback filter (FBF) but any residual uncausal part will increase the MSE, causing 
performance degradation. Enlarging the value of N will decrease the residual uncausal part 
at the price of increased computational complexity. 

In the embodiment of FIGURE 4, the C coefficient calculator 41 finds the WLS 
solution for the number of coefficients chosen to be initialized, using either Equation 4 above 
for a T-spaced equalizer or Equation 5 above for a fractionally spaced equalizer, and using 
the exemplary linear system and the exemplary weight matrix described below. 

The equations of the linear system in this example are: 

Equation 6: 

V 

Z dk]H[-k]=l 
k=-U 
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Equation 7: 

Z^C{k]-H{n- k] = 0, #i = -1-2-3... 

Equation 8: 

1 C\k] H[n- k] = 0, w = -2-4... 

*=-t/ 

Equation 9: 

E C[Ar] = 0, w = 12,3.- 

Equation 10: 

Z <M •//[/*-£]=(),« =2,4,.. 

Equation 1 1 : 
C[£]=0>k=l,2,3...K 

The foregoing equations of this overdetermined linear system (i.e., having more 
equations than variables) basically represent the following conditions. Equation 6 reflects 
that the lead tap is 1 in the equivalent response of H[n] and C[n]. Equations 7 and 8 relate 
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to uncausal ISI cancellation, and basically reflect zero uncausal ISI in the equivalent response 
of H[n] and C[n], Equation 7 corresponding to the T-spaced equalizer and Equation 8 
corresponding to the fractionally spaced equalizer. Equations 9 and 10 relate to causal ISI 
cancellation, and basically reflect zero causal ISI in the equivalent response of H[n] and C[n], 
Equation 9 corresponding to a T-spaced equalizer and Equation 10 corresponding to a 
fractionally spaced equalizer. Equation 1 1 basically reflects minimum noise gain by C[n]. 

An exemplary weight matrix for use in the WLS solution will now be described. 
Each of Equations 6- 1 1 actually represent one or more equations. Referencing Equation 6, 
there is clearly only one lead-tap equation. With reference to the uncausal ISI cancellation 
10 j (S & addressed by Equation 7 (or Equation 8), let the number of equations actually represented 
□ by Equation 7 (or 8) be designated as p. Referencing the causal ISI cancellation addressed 
j ;3 by Equation 9 (or 10), let the number of equations actually represented be designated as q. 
The p and q values are generally implementation dependent and can be set as desired. In 
practice, the designer considers the number of taps (causal/uncausal) of the H[n]*C[n] 
15 response, where * is the convolution operator. If the number of uncausal taps in H[n] is T 

and the number of uncausal taps in C[n] is U, then the number of uncausal taps in H[n]*C[n] 
is U+T or less. Regarding the minimum noise gain addressed by Equation 11, there are K 
equations because, in this example, C[n] has K coefficients. 
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With the aforementioned definitions of p, q and K, the weight matrix is given by: 
Equation 12: 



W = 



a L 



where \ denotes the unit matrix of size j x j. If the noise variance, NO/2 is given (if it is not 
given it can be conventionally estimated), the minimum MSE can be achieved by setting the 
following relationships: 
Equation 13: 

a = 0 (or q = o) 



10 



Equation 14: 



b = 



2 < 
NO 



Equation 13 means that the causal ISI cancellation equations (corresponding to Equations 
9 or 10) will not be used. Also note that the weight b of Equation 14 is the reciprocal of the 
signal-to-noise ratio (SNR). 
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In some embodiments, it is important to reduce the size (i.e., the absolute value) of 
the D coefficients, at the price of increased MSE. In such embodiments, the causal ISI 
cancellation equations should be used, thus requiring a non-zero value of the weight a. For 
example, a = 1 gives causal ISI cancellation the same weight as uncausal ISI cancellation. 
Other examples include a = 0.1 and a = 0.5. 

Using Equations 4-14 appropriately for either a T-spaced equalizer or a fractionally 
spaced equalizer, the C coefficient calculator 41 can find the WLS solution for the selected 
number of coefficients of C which are to be initialized. Calculation of the WLS solution for 
the selected number of coefficients using Equations 4- 1 4 is a straightforward process that can 
be readily implemented by workers in the art, as described in detail in H. L. Van-Trees, 
Detection. Estimation and Modulation Theory, Part 1 , (Wiley, 1968), and H. W. Sorenson, 
Parameter Estimation , (Marcel-Dekker, 1980), both of which are incorporated herein by 
reference. 

The initial subset of C coefficient values produced by the C coefficient calculator 41 
is provided (together with the aforementioned channel estimation information) to the D 
coefficient calculator 43, which calculates in response thereto the D coefficient values. The 
D coefficient values calculated by the D coefficient calculator 43, together with the subset 
of C coefficient values calculated by the C coefficient calculator 41 , are input to the selector 
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37 of FIGURE 3, which in turn provides this initial subset of the equalizer filter coefficients 
to the feed forward equalizer 13 (initial subset of C coefficients) and the feedback filter 17 
(initial D coefficients). After the D coefficient calculator 43 calculates the initial D 
coefficient values and passes them through the selector 37 (along with the initial subset of 
C coefficient values), the control signal 36 switches the selector so the subsequent 
coefficients for C and D are thereafter provided from output 38 of the coefficient adaptor 35 
of FIGURE 3. 

In response to the initial subset of coefficient values of C produced by the C 
coefficient calculator 41 , the exemplary D coefficient calculator 43 in FIGURE 4 calculates 
the initial D coefficient values (a total of L values) using the following exemplary equation 
(for a fractionally spaced equalizer). 

Equation 15: 



FIGURE 5 illustrates exemplary operations which can be performed by the 
embodiments of FIGURES 3 and 4. When it is determined at 50 that a communication burst 
has arrived, the taps of H[n] (i.e., lead, causal and uncausal) are assigned at 52. At 54, initial 
values of a subset of the coefficients of C[n] are calculated according to a first algorithm 
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(e.g., LS or WLS). Thereafter at 56, the initialized coefficient values of C[n] are used along 
with the known H[n] information to calculate initial values of the coefficients of D[n]. 
Thereafter at 58, subsequent values of all the C[n] and D[n] coefficients are calculated based 
on the initial values, using a second algorithm (e.g., LMS). The subsequent coefficient 
values are then used at 59 to decode the unknown symbols. The operations at 58 and 59 are 
repeated for every M unknown symbols, for example, for every unknown symbol (M=l), 
until the end of the burst at 60. 

It can be seen from the foregoing description that the present invention utilizes a first 
coefficient determination technique to provide an initial estimate of the values of a subset of 
the feed forward coefficients and, from this subset of feed forward coefficient values, 
calculates initial values of the feedback coefficients. This initialization of the feedback 
coefficients and a subset of the feed forward coefficients provides an initial slicer MSE that 
is low enough to permitted subsequent symbol decoding using a second coefficient 
calculation technique that is less computationally complex than the first coefficient 
calculation technique. This can advantageously reduce both latency, and the number of 
MMACs required for equalizer operation, while also providing an initial slicer MSE that is 
sufficiently low to permit decoding with a desired symbol error rate. 
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It will also be evident to workers in the art that the embodiments of FIGURES 3-5 
can be readily implemented, for example, by suitable modification of software, hardware, or 
a combination of both, in the equalizer portion of a conventional digital communication 
receiver. 

Although exemplary embodiments of the invention are described above in detail, this 
does not limit the scope of the invention, which can be practiced in a variety of 
embodiments. 
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